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ABSTRACT 

The  dynamic  simulation  of  multibody  tracked  vehicles  offers  engineers  a  powerful  tool 
with  which  they  may  analyze  and  design.  Currently,  parts  of  these  complex  mechanisms  are 
introduced  to  multibody  algorithms  as  rigid  bodies.  Then  in  a  follow-on  structural  analysis,  the 
loads  from  the  multibody  dynamic  simulation  are  input  to  calculate  strains  and  stresses  within  the 
bodies.  The  present  investigation  aims  to  establish  appropriate  means  by  which  flexible  three- 
dimensional  track  links,  which  allow  large  relative  rotation  between  the  elements,  can  be 
modeled.  This  will  pave  the  way  towards  the  incorporation  of  detailed  flexible  structural  models 
into  a  multibody  dynamic  simulation  environment  allowing  for  an  integrated  solution.  In  addition, 
a  new  formulation  for  the  interaction  between  the  rigid  sprocket  teeth  and  flexible  chain  is 
presented.  Numerical  results  are  introduced  to  illustrate  the  effects  of  flexible  links  on  the 
dynamics  of  tracked  vehicles. 


INTRODUCTION 

Current  commercial  implementations  for  tracked  vehicle 
modeling  are  not  ideal  for  the  analysis  of  three-dimensional 
flexible  chain  tracked  vehicles.  Each  track  chain  may  contain 
large  number  of  separate  bodies/elements  which  require 
connectivity  conditions  be  imposed  at  each  time  step. 

Non-multibody  simulation  models  have  had  success  in 
providing  engineers  with  a  means  with  which  to 
quantitatively  compare  tracked  vehicles  [1-3].  It  is  hoped, 
however,  that  multibody  simulation  models  will  provide  the 
working  engineer  tracked  vehicle  models  based  on  fewer 
assumptions  and  greater  flexibility  of  simulation. 

It  is  one  of  the  objectives  of  this  investigation  to  present 
the  results  of  incorporating  the  recently  developed  Absolute 
Nodal  Coordinate  Formulation  (ANCF)  C'VC1  flexible  chain 
into  a  tracked  vehicle  multibody  model.  The  C°/C’  chain 
model  offers  a  number  of  advantages  over  current  chain 
models.  First,  the  linear  connectivity  conditions  resulting 
from  ANCF  may  be  applied  at  a  preprocessing  stage  [5],  This 
allows  for  an  efficient  elimination  of  the  dependent  variables. 
When  used  in  conjunction  with  ANCF  spatial  finite  elements, 
the  linear  connectivity  conditions  lead  to  a  constant  inertia 


matrix,  as  well  as  zero  Coriolis  and  centrifugal  forces. 
Further,  instead  of  considering  a  series  of  interconnected 
rigid  or  flexible  bodies/elements,  the  new  formulation  allows 
for  the  use  of  one  finite  element  mesh  for  the  chains  despite 
the  finite  rotations  of  the  track  links  with  respect  to  each 
other  and  relative  translation  between  the  two  track  chains. 
Numerical  examples  will  be  used  to  illustrate  the  advantages 
of  the  C’VC1  flexible  chain  formulation  as  relevant  to  tracked 
vehicle  models. 

Another  objective  of  this  investigation  is  to  present  a 
newly  developed  formulation  for  the  contact  between  the 
C'VC1  chain  model  and  the  sprocket  teeth.  This  interaction 
formulation  will  allow  the  building  of  tracked  vehicle  models 
that  can  more  accurately  represent  the  complex  dynamics  of 
tracked  vehicles. 

The  organization  of  the  paper  is  as  follows.  The  first 
section  contains  a  summary  of  the  absolute  nodal  coordinate 
formulation  (ANCF),  as  well  as  a  presentation  of  the  beam 
element  used  in  this  investigation.  In  the  next  section,  the 
constrained  multibody  dynamic  equations  of  motion  are 
presented.  The  third  section  summarizes  the  recently 
developed  finite  element  ANCF  mesh  which  can  be  used  to 
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model  flexible  chains.  The  fourth  section  presents  the 
models  used  for  the  interaction  between  the  different  parts 
of  the  tracked  vehicle  running  gear.  Specifically,  the  fourth 
section  will  introduce  a  new  formulation  for  the  simulation 
of  sprocket-flexible  chain  interaction. 


1.  ABSOLUTE  NODAL  COORDINATE 
FORMULATION 


In  this  investigation,  the  finite  element  absolute  nodal 
coordinate  formulation  is  used  in  the  modeling  of  flexible 
bodies.  The  global  position  of  an  arbitrary  point  on  element 
e  of  body  i  can  be  defined  (figure  1)  in  an  inertial  global 
coordinate  system  XYZ  as 

r  =S  (x  ,y  ,z  )e 

where  S,e(xie ,  yle ,  ze)  is  the  space-dependent  matrix  that 

defines  the  element  shape  functions,  and  e,cis  the  time- 
dependent  vector  of  nodal  coordinates. 

The  vector  of  nodal  coordinates  e“  consists  of  nodal 
displacements  and  gradients.  It  can  be  defined  as 


e"=[r"( 


c)r'H  })rA  r)r" 


r)r"2  c)r"2  drk 


r)x'L  df  c)z,c  dx‘e  df  r)zK 


-f(2) 


Beam  centerline 


Figure  1:  ANCF  Beam  Element  Coordinates 

In  equation  (4),  £  =  Je / / ,  tj  =  y*e// ,  g  =  z*e/ / ,  and  1  is 
the  length  of  the  element.  It  can  be  shown  that  the  three 
dimensional  beam  element  presented  here  is  capable  of 
capturing  the  effects  of  shear  deformation  and  rotary  inertia 
as  well  as  the  coupling  between  the  cross  section 
deformation  and  bending/extension. 


„iek  1  iek  1  M  iek  1 

In  equation  (2),  riekT  - —  - — ,  and  - — —  refer 

a  x  d  y  d  z 

to  the  displacements  and  gradients  at  node  k.  It  can  be  shown 
that  the  selection  of  nodal  positions  and  gradients  as  element 
coordinates  leads  to  an  exact  modeling  of  rigid  body 
dynamics.  The  shape  function  for  the  three-dimensional 
beam  element  based  on  the  absolute  nodal  coordinate 
formulation  used  in  this  investigation  is  defined  by 

S'C  =  [.sfl  sl I  si  I  si  I  sfl  si  I  si  I  s‘eI]  (3) 

where  I  is  the  3x3  identity  matrix  and  the  shape  functions 
s'i ,  k  =  are  given  as  follows  [4]: 

5,=l-3f +  2f,  S2  =  ltf- +  ^3),  S3  =  im  1-4) 

S4  =  Ig (1  -  a  S5  =  3£3  -  2f,  S6  =  l(-f  +  ^3)  (4) 

S7  =  Ss  =  lgZ 


2.  CONSTRAINED  EQUATIONS  OF  MOTION 

Using  the  principle  of  virtual  work,  one  may  write  the 
finite  element  equations  of  motion  in  terms  of  the  vector  of 
nodal  coordinates  as  follows: 

Mie&ie  =(£ +Ql  (5) 

In  equation  (5),  Q)  is  the  vector  of  generalized  element 

elastic  forces,  Q“  is  the  vector  of  generalized  element 

external  forces,  and  M,e  is  the  constant  mass  matrix  of  the 
element. 

Using  equation  (5),  one  may  write  the  equations  of  motion 
of  the  constrained  multibody  system  in  the  following  form: 

C(q,f)  =  0  j 
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where  M  is  the  system  mass  matrix,  q  is  the  vector  of  the 
system  coordinates  (which  includes  the  absolute  nodal 
coordinates  and  other  coordinates  used  to  describe  the 
motion  of  the  rigid  bodies),  X  is  the  vector  of  Lagrange 
multipliers,  Q  is  force  vector,  C  is  the  vector  of  kinematic 

constraint  equations,  and  Cq  is  the  Jacobian  matrix  of 
constraint  equations. 


3.  NEW  ANCF  FINITE  ELEMENT  MESH 

Fully  parameterized  absolute  nodal  coordinate  formulation 
finite  elements  have  12  degrees  of  freedom.  These  degrees 
of  freedom  include:  three  rigid  body  translations,  three 
rotations,  and  six  deformation  modes.  We  may  introduce  a 
new  formulation  for  a  pin  joint  that  uses  the  following  six 
scalar  equations  defined  at  the  joint  node  between  bodies  i 
and  /  [51: 

r  =  rJ  (7) 

3r'  _  drj  (8) 

da  da 

In  the  above  equations,  a  is  the  coordinate  line  that  defines 
the  joint  axis.  It  can  be  either  x,  y,  z,  or  any  other  coordinate 
line  since  the  use  of  ANCF  makes  easy  the  use  of  the 
gradient  tensor  transformation.  Consider,  for  example,  the 
gradient  transformation  between  u,  v,  w  and  x,  y,  z.  This 
transformation  is  given  by  [r„  rr  rH,]  =  [rr  rv  r  Ja 
where  A  is  given  by  the  constant  matrix  shown  below. 


dx 

dx 

dx 

du 

dv 

3  w 

dy_ 

du 

3v 

3  w 

dz 

3z 

dz 

du 

3v 

3  w 

The  six  scalar  equations  represented  by  equations  (7)  and 
(8)  reduce  the  number  of  degrees  of  freedom  for  the  new  pin 
joint  to  six;  one  rotation,  and  five  deformation  mode  degrees 
of  freedom.  As  a  result,  we  have  C1  continuity  with  respect 
to  the  a  coordinate  and  C°  continuity  with  respect  to  the 
remaining  coordinates.  The  linear  scalar  equations  (7)  and 
(8)  can  be  imposed  at  a  preprocessing  stage,  thus  allowing 
for  the  elimination  of  the  dependent  variables.  These 
equations  can  be  used  to  develop  a  new  kinematically  linear 
finite  element  mesh.  This  would  allow  arbitrarily  large 


relative  rotations  between  its  constituent  elements  at  the 
mesh  joints. 


Figure  2:  Geometry  of  contact  between  rigid  sprocket 
teeth  and  the  C°/C 1  chain 


4.  FLEXIBLE  CHAIN-RIGID  COMPONENT 
CONTACT  FORMULATIONS 

The  new  finite  element  ANCF  C'VC 1  mesh  formulation  that 
was  presented  in  the  previous  section  can  be  used  to  more 
accurately  describe,  with  one  mesh,  a  flexible  track.  In  this 
section  the  contact  between  the  flexible  track  links, 
sprockets,  rollers,  and  idlers  will  be  described.  The  idlers, 
sprockets,  and  rollers  will  be  modeled  as  rigid  bodies  in  this 
investigation. 

While  previous  work  has  focused  on  modeling  flexible  belt 
driven  vehicles  [6],  the  present  investigation  differs  in  that  it 
includes  a  formulation  for  the  interaction  between  the 
flexible  chain  and  sprocket  teeth.  This  will  allow  for  future 
investigations  in  which  the  dynamics  of  rigid  and  flexible 
link  tracked  vehicles  can  appropriately  be  compared. 

The  formulation  for  the  interaction  between  the  sprocket 
teeth  and  the  C°/C'  chain  is  as  follows  (one  can  refer  to 
figure  2  for  the  definitions  of  the  variables  used):  let  N  be 
the  vector  that  defines  the  local  position  of  the  center  of  the 
pin  joint  in  the  sprocket  coordinate  system,  U  is  a  vector 
which  defines  the  axis  around  which  the  sprocket  rotates, 
r  ■  is  the  pin  radius,  lg  is  the  backlash  as  a  percentage  of 

the  pin  radius  r  ■  ,  n,  is  the  current  tooth  number,  and  T  is 
the  number  of  teeth.  Note  that  it  is  assumed  that  the  sprocket 
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axis  is  parallel  to  the  Y-axis.  Take  angle  a  as  from  the  center 
line  of  the  tooth  to  the  vertical  center  line  of  the  sprocket. 
The  two  angles  a  and  p  can  be  calculated  from  geometry  as 

/  \  2  71 

a  =  (nt  —  l) —  (10) 

T 


P  =  tan 


-1 


( l  +  lB ) 


pin 


N 


(11) 


With  the  tooth  contact  force  established,  one  also  needs  to 
include  the  contact  force  due  to  the  sprocket  base  circle  and 
pin  contact.  Figure  3  illustrates  the  geometry  used  to  define 
the  contact  between  the  C°/C'  chain  and  the  cylindrical  rigid 
body  that  constitutes  the  sprocket  base  circle  [6]. 

The  radial  (normal)  penetration  for  the  configuration 
shown  in  figure  3  can  be  formulated  as  follows: 


One  can  then  find  the  vectors  that  define  the  location  of  the 
contact  points  between  the  sprocket  teeth  and  the  pin  as 


pi 

vp2- 


fpl 

p2 


[cos(a+/?)  0  sin(a+ ft)  ]^ 

j^cos  {a-p)  0  sin  {a-p)]  T  (13) 


where  the  magnitude  of  r^andr^  can  be  calculated  from 
equation  (14). 

>1 1  =  |rp 2  |  =  ^N2+(l  +  /B)2  (V„)2]  (14) 

The  distance  between  the  pin  and  the  tooth  contact  point  is 
given  by  equations  (15-16). 


Vi  =  Vi  -  N 

rn2=rp2~^ 


The  condition  for  contact  is  given  by 


(15) 

(16) 


r  <  r  or  r  „  <  r . 

ln  1  —  ,pin  WA  ln2  —  'pin 


(17) 


Figure  3:  Geometry  of  contact  between  rigid  cylinder 

(21) 


and  the  C°/C'  chain 


d  =  Rc~  |N| 

where  R  is  the  radius  of  the  cylinder  and  N  is  a  vector  that 
defines  the  position  of  the  contact  point  in  the  coordinate 
system  of  the  cylinder.  N  can  be  defined  as 


N  =  r  - 


(r;fc)' 


(22) 


If  the  conditions  for  contact  are  true,  the  force  normal  to 
the  contact  surface  is  evaluated  using  a  continuous  force 
model  defined  as  follows 

F n=(kd  +  cd)V  (18> 

where  d  can  be  found  from  the  following  equation 

d  =  rpin-\rnl\  or  d=  |r«2 1 —  rpin  <19> 


The  normal  force  for  the  contact  between  the  rigid  cylinder 
and  the  C^C1  chain  can  be  then  be  found  by  using  the 
appropriate  penetration  and  rate  of  penetration  in  equation 
(18). 

The  tangential  contact  force  is  based  on  a  trilinear 
Coulomb  friction  law  characterized  by  the  following 
equations 

Ft  =-//|F„|t,  |  Vt  |  >  v0  (23) 

Ft=-vJVt|t,  |Vt|<v0  (24) 


and  the  vector  normal  to  the  contact  surface  is  easily  found 
from 


V  =  U  x  N 


(20) 


where  V  is  the  tangential  velocity,  V  is  a  friction 
parameter  that  defines  the  slope  in  the  transition  area,  and 
the  constant  V0  is  defined  as 


V 


o 


(25) 
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The  total  force  will  then  be  the  summation  of  the  normal 
and  tangential  forces.  In  the  case  of  the  idlers  and  rollers,  the 
contact  forces  will  be  calculated  in  a  manner  similar  to  the 
contact  between  the  sprocket  base  circle  and  the  C'VC1  chain. 


5.  NUMERICAL  RESULTS 


The  previous  section  presented  the  formulations  used  in 
modeling  the  contact  between  the  rigid  sprockets,  rollers, 
idlers,  and  the  C'VC1  track  or  chain.  In  this  section, 
numerical  results  of  a  model  of  a  rubber  tracked  vehicle  that 
includes  only  one  track  are  presented.  The  model  contains  a 
ground,  chassis,  sprocket,  idler,  two  rollers,  and  a  flexible 
track.  The  mass  of  the  driving  sprocket  and  the  idler  are 
assumed  to  be  40  kg.  The  inertia  of  driving  sprocket  and  the 
idler  about  the  X,  Y,  and  Z  axes  are  respectively  0.4  kg.nr2, 
0.8  kg.nr  and  0.4  kg.nr2.  The  mass  of  the  rollers  is  10  kg. 
The  mass  of  the  vehicle  chassis  is  100  kg.  In  the  numerical 
example  presented  in  this  section,  the  chassis  is  assumed  to 
have  zero  initial  velocity.  The  angular  velocity  constraint 
that  is  applied  to  the  sprocket  is  defined  by  the  following 
equation: 


Co 


t  <=0.2 


&>=  < 


100 


Q-0.2) 

0.8 


100 


0.2  <t<=  1  (26) 

t  >0.2 


7  ,  3  , 

respectively,  by  9x10  N/nr  and  2x10  N.s/m  .  Tangential 
friction  forces  are  also  introduced  using  a  coefficient  of  dry 
friction  of  1.2.  The  friction  parameter  that  defines  the  slope 
in  the  transition  region  is  assumed  to  be  106  N.s/m3.  Figure 
4  shows  the  X,  Y  and  Z  coordinates  of  the  position  vector  of 
the  center  of  mass  of  the  chassis.  Figure  5  shows  the  velocity 
of  the  chassis  in  the  X-direction.  Note  that  in  figure  6  the 
displacement  decreases  at  certain  times;  this  is  an  indication 
of  sliding  behavior,  which  will  continue  until  the  rubber 
elements  in  contact  with  the  ground  obtain  zero  speed  in  the 
x  direction. 


where  cn  is  the  rotational  speed  of  sprocket  in  rad/s  and  t  is 
the  time  in  seconds. 

The  diameter  of  the  base  circle  of  the  sprocket  is  0.38  m 
and  has  10  teeth.  The  diameter  of  the  idler  and  the  rollers  are 
respectively  0.38  m  and  0.18  m.  The  chain  is  modeled  using 
24  ANCF  elements  that  have  a  rectangular  cross  section  of 
dimensions  0.02mx0.4m  and  density  of  2000  kg/m3.  The 
incompressible  Neo-Hookean  model  and  non-linear  damping 
model  are  used  to  model  the  internal  behavior  of  the  rubber 

g 

chain.  The  Neo-Hookean  constant  is  assumed  to  be  1x10  , 

9 

the  incompressibility  constant  is  assumed  to  be  1x10  ,  and 
the  dilatation  and  deviatoric  dissipation  factors  used  for  the 

-4  -5 

damping  model  are  assumed  to  be  1x10  and  5x10  , 

respectively.  A  pretension  positive  axial  strain  of  value  0.05 
is  applied  to  the  rubber  chain.  A  compliant  force  model  is 
used  to  describe  the  contact  between  rigid  bodies  and  the 
chain.  The  stiffness  and  damping  coefficients  used  in  the 
rigid  bodies  and  chain  contact  force  model  are  given, 


Figure  4:  Position  of  Center  of  Mass  of  the  Chassis 
( — • —  X  coordinate;  — ■ — Y  coordinate;  ^  Z 

coordinate) 


Figure  5:  X  velocity  of  the  center  of  mass  of  the  chassis 
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Figures  6  and  7  show,  respectively,  the  X  and  Z 
coordinates  of  a  point  on  the  rubber  chain  that  is  initially 
located  at  the  top  of  the  idler,  as  function  of  time.  Figure  8 
shows  the  motion  trajectory  of  this  point  relative  to  the 
moving  center  of  mass  of  the  chassis. 


Figure  6:  X  coordinate  of  the  position  vector  of  a  point 
on  the  rubber  chain 


Figure  7:  Z  coordinate  of  the  position  vector  of  a 
point  on  the  rubber  chain 


X-position  (m) 

Figure  8:  Trajectory  of  the  position  vector  of  a  point  on 
the  rubber  chain 

6.  SUMMARY  AND  CONCLUSIONS 

The  present  investigation  proposed  a  method  to  model  the 
contact  between  spatial  track  chain  links  and  rigid  running 
gear  components.  The  Absolute  nodal  coordinate 
formulation  was  used  to  model  the  C°/C’  track  or  chain 
using  linear  connectivity  conditions.  This  chain  formulation 
allows  for  C1  continuity  in  the  direction  of  the  pin  joint 
definition  and  C°  continuity  with  respect  to  the  other 
components.  A  new  formulation  for  the  compliant  force 
contact  between  three  dimensional  flexible  track  links  and 
rigid  sprocket  was  presented.  Along  with  the  formulation  for 
the  contact  between  the  rigid  rollers,  idlers,  and  flexible 
links,  these  contact  formulations  were  used  to  model  a 
tracked  vehicle  with  a  single  track.  This  model  may  then  be 
used  in  future  investigations  to  compare  the  dynamics  of 
rigid  and  flexible  tracks  in  tracked  vehicle  models. 
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